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1 . Abstract 


We have extended the recent work of Allen, Gayle, and 
Richardson (1970) and Suits (1972) to compute directional 
reflectance from a row crop. We assume a 'row canopy' to he 
bounded in the horizontal direction perpendicular to the rows, and 
in the vertical direction while along the rows the 'row 
canopy' is assumed to extend to infinity. The problem can 
be easily adapted to a crop in which the row structure is 
destroyed by the components of plants in one row running 
into components of another row. In this case one assumes 
that the plant canopy extends to infinity in the horizontal 
direction perpendicular to the rows. The advantage of 
bounding a plant canopy in horizontal directions is that 
the edge effects are properly taken into account. The 
calculations-'can easily be extended to incorporate a boundary 
between two fields. The incident solar radiation is de- 
composed into three cartesian components, one along the 
rows and two perpendicular to row direction. For horizontal 
incidence the scattering and absorption coefficients of 
specular radiation are assumed to be different from those for 
vertical incidence. Hopefully, the introduction of these new 
scattering coefficients will help in discriminating crops with 
leaf orientation in preferred directions. Changes in leaf 
orientation caused by certain diseases will also show up as 
changes in reflectance and transmittance. 



This report only contains development of mathematical 
equations. Numerical results based upon these equations will 
be published in a forthcoming report. 



Introduction . 


Attempts at identifying vegetation types have been 
made in more than one way. Laboratory experiments per- 
formed to study the absorption spectra of a single leaf 
taken from each vegetation type(see Me Leod, 1971) have 
not yielded algorithms suitable for discriminating one 
type of vegetation from another. The main reason for 
this is that leaves from most plants contain almost the 
same amount and type of chlorophyls and wa ter ,~ which- are 
responsible for the absorption spectra. Another approach 
that has been applied to this problem is to study the scat- 
tering and absorption properties of a collection of fresh 
leaves placed flat on a background of known reflectivity. This 
approach represents an improvement from the one leaf case as 
this arrangement is closer to a plant than an individual leaf. 
However, it is still not close enough to a vegetation field 
or even to a plant. 

The aim in a mathematical model of a plant canopy is the 
same as in the laboratory experiments, that is to arrive at 
algorithms which will help in discriminating vegetation 
types. The mathematical models developed in the past assume 
a plant to be occupied between two horizontal planes, extending 
to infinity in all horizontal directions. These models 
can only give reliable results if the sensors view area is 



small enough, that is, either the plant canopy is large 
or a vegetative field being observed has closely spaced 
plants and a small portion of the field is being viewed 
at a time. Furthermore, if the field of view contains 
more than one vegetation type the canopy model needs mod- 
ification in that boundaries in horizontal directions 
need t'o be added. 

When the sun is at an angle from the local vertical 
the light penetrates through a horizontal direction in 
addition to the vertical direction. Now, consider that 
the object under observation is a row crop and the sun 
is at a non zero angle from the local vertical and the 
azimuthal angle of the sun with respect to the row direc- 
tion is also non zero. It is obvious that for a complete 
description of the radiation field we need to solve for 
the radiation equations in three mutually orthogonal 
directions, one vertical and two horizontal directions 
parallel and perpendicular to the rows. Section 3 of this 
paper is devoted to this problem. In sections 4, 5, and. 6 
we shall elaborate on additional shadowing of soil between 
plants, limitations of our model and some experiments that 
should be performed to gain insight into, plant canopy ef- 
fects. This report only contains development of mathe- 
matical equations . ^ Numerical results based upon these 

equations will be published in a forthcoming report. 



3 . Review of Previous Plant Canopy Models 

Allen and Richardson (1968) were the first to employ 
the equations of Kubelka and Munk (referred to below as KM) 
(1931) to a plant canopy. The KM and a large number of 
related papers were concerned with the study of interaction 
of light with paints, glass, paper, and plastic materials 
where the object consists of tightly packed particles. In 
general light penetrates only a short distance into.the 
material and the lateral extent of the material does not 
play an important role. For this reason it sufficed to 
assume that the object extended to infinity in lateral 
direc tions . 

The only scattering coefficients occuring in the KM 
work are an absorption coefficient and a back scattering 
coefficient, both of these coefficients are related to in- 
ternal diffuse light. Duntley (1942) reported theoretical 
work of Ryde (1931,1932) who, based on comparison of his 
work to observed data, concluded that absorption and scat- 
tering coefficients of incoming specular light are different 
from those of internal diffuse light(Since these equations are 
now being used in plant canopy models it seems appropriate to 
verify this same result experimentally for plant canopies). The 
mathematical equations of Duntley (1942) which incorporate 
this modification contain five unknown constants, two for 
diffuse light and three for specular. Allen, Gayle and 
Richardson (1970) used these improved equations to study 



plant canopy effects. Again, as with the KM work, the 
plant canopy is assumed to have infinite lateral extent. 

Two important assumptions in the KM work are: (1) Lambert 

Cosine law, i.e. there is no Fresnel reflection. It is not 

known to the present author if this assumption has been 
experimentally verified in the case of plants and leaves. 

(ii) the particles in the layer are regarded as randomly 
distributed and smaller in size than the thickness of the 
layer. Referring to the first assumption it is quite pos- 
sible that some percentage of the plant area acts as a spec- 
ular reflector. However, because of lack of time this sub- 
ject will not be dealt with in this report. Concerning 
the second assumption, one observation can be made by visual 
examination, that is, in a healthy plant the. upper surface 
of a leaf is, in general, a better reflector than its lower 
surface. This fact will be incorporated in the present work. 
Suits (1972) made a useful extension of the work of Allen, 
Gayle and Richardson in that he defined a vertical leaf 
area index of a plant. The leaves were assumed to be 
Lambertian reflectors but the introduction of the new ver- 
tical leaf area index leads to expressions of transmitted 
and reflected radiation which depend upon the sun and sen- 
sor angles, thereby yielding non-Lamb e r t.ian expressions for the 
reflectance and transmittance. However, Suits has assumed 
that the reflection and transmission coefficients for. 
specular and diffuse radiation are identical, as can be seen 
from his expressions (6) and (7). 



4. A Model for Interaction of Light with a Row Crop 


In this section we shall obtain a self consistent 

radiation field in a row crop whose plant height, row width 

and row direction are known. Also known are the angular 

coordinates (9 , i|r ) of the sun where 9 represents the zenith 

s s s 

angle of the sun and ijr g its azimuth with respect to the 

row direction. If Iq is the intensity of sunlight then 

the components along the downward vertical (z-axis), along 

the rows (xj-axis) and pe rp ind icu 1 ar to the rows (y-axis) 

are I n Cos9 , I n Sin9 Cos<|f and IS in9 Sin>|> respectively. 

The coordinates and components are shown in figure 1. 

Let us denote the attenuated components of specular light 

as it passes through the canopy by I , I , I and the 

s x s y s z 

components of diffuse light in the positive and negative 

directions along the coordinate axes as I, , I , I, , I , 

+x - x +y - y 

I , and I z respectively. In the' vertical direction we 

introduce |A z> B^, as absorption, back scattering and 

forward scattering coefficients for specular light. The 

quantities u , B , F and la , B , F are similar coefficients 
M x x x y y y 

for.the x and y components of specular light. Most plants 

have axial symmetry and for this reason we shall assume 

la =u , B =B , F =F . For diffuse light in the z direction 
r x r y’ x y x y 

we assume back scattering coefficient 'B' of I + ^ to be 
different from ' 1$ ' the back scattering coefficient for I ^ ; 
this we introduce to account for greater reflectance from the 
upper surfaces of leaves. For the x and y components of 




diffuse radiation the back scattering coefficient is assumed 
to be B" different from.B or B’ as the scattering in horizontal 
directions depends upon the vertical ^leaf area index which can 
be different in magnitude from the horizontal leaf area index. 
Similarly the absorption coefficient for x,y component 
equations is assumed to be a new parameter ^ ' . With these 
approximations the equations for specular and diffuse light 
in the three coordinate axes directions can be written as: 


Z - Ax i s 


d I 

s z 
dz 


dl 


+z 


dz 


- (M> +B +F ) I 
z z z s z 


FI - ( p, +B ) I + B 

z sz +z 


I 

- z 


( 1 ) 

( 2 ) 


dl 

- z 
dz 


B I 
z s z 


(p+B ' ) I + BI, 

- Z +Z 


( 3 ) 


X- Axis : 


dl 


sx 


dx 


- (P +B +F )I 

X X X sx 


(A) 


dI +x 

3 = F I - ( p+B " ) I + B " I 

dx x sx +x -x 


( 5 ) 


dl 

- x 
dx 


B I 

X S X 


(p+B " ) I + B " I 

- x +x 


Along the y-axis equations can be written simply by 
replacing I gx by I sy , P x by p- y 


( 6 ) 


etc. The boundary con- 



ditions along the x-,y-,z-axes are all different. For ex- 
ample, at the top of a plant, which is at z = 0 and at the sun- 
lit side of the row i.e. y=0, and at the sunlit end of a row 
there is no diffuse light in the positive directions; i.e. 

(7) 

( 8 ) 

(9) 

At the bottom of a row we have 

I - z ■ R *< 1 +, M „) at z ' z i ■ < 10 > 

At the other side of the row 

I = 0 at y=y 1 (11) 

where y is the width of plant cover in a row. At the other 
end of a row (which is assumed to be of infinite length) 

I_ =0 at x=°° (12) 

The solution of equation (l)-(6) and the y-equations sub- 
ject to boundary conditions (7)- (12) need not be written 
down as we are only interested in quantities that can be 
observed with the help of a remote sensor. These quanti- 
ties are the reflectance and transmittance at five surfaces: 
x = 0 , x=°o, y = 0 , y=y^, and z = 0 . Reflectance values at z = 0 , x = 0 

and y=0 are given by I (z=0), I (x=0) and I (y = o) 

“ z “X - y 

respectively. Transmittance at x=<* ig 0 and fchdt'at y =y ^ is 


l +z (z=0) = 0 

i +y (y=o) = o 

l +x (x=0) = o 
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I + y(y=y^). These quantities are given below: 


I_ 

(2 = 0) 
z 

- A 3 + A 4 + N l iQ 

(13) 

whe re 

X z 0 

= I 0 COS 8 S 

(14) 


Q 3 = 

M- + B + F 
Z Z Z 

(15) 


A 3 " 

^ +B+a l A = 
B ' 

(16) 


A 4 = 

“ +B+ “2 A =MA, 
B ' 

(17) 



a 2 Z l 

PI z0 (M-Rg)e z + (l-PRg-N)I zl 

(18) 


A i 

a i z i a 2 Z l 

e 1 1 (L-Rg) - (M-Rg) e z 1 


X 21 = 

1 z0 «=’'P<-Q 3 ^ 1 ) 

(19) 


A = 

a z 

-PI z0 (L-Rg)e 1 1 - (l-PRg-N)I zl 

(20) 


A 2 

N = 

CC]Z. a 2 Z T 

e 1 1 (L-Rg) - (M-Rg) e z 1 

(p-+B-Q Q )B +BF 
J z z 

(21) 


(li-Oj) (|i-a? 2 ) 


P = 

B'B + (p.+B '+Q,)F Z 
z J 

(22) 


C 

(^-0^) . (M.-Qf 2 ) 


“l = 

-i(B-B') + Jp, 2 +y,(B+B ' )+£(B-B ' ) 2 

(23) 


“2 = 

-|(B-B') - Jm- 2 +P.(B+B ' )+£(B-B ' ) 2 

(24) 

and Rg 

is the 

reflectivity of the ground. 



i (y 
+y 

. . Py i " Py l A NUM 

= y ! ) = A 12 e + A 22 e + Den 3 

- 1 y i < 25 > 


whe re 



P 2 X V Q M 2 e 1 + 

a 12 py ■, -Py 


l rj l 

L 2 e - M 2 e 


2’ y 1 s ANUM 2 
DEN 4 


( 26 ) 


8y 

'P,I.nL, e “ (l-No)I 


22 


2 yO^ 


2 y 1 


DEN 4 


(27) 


B "B 2 + (M,'+B"+Q 2 )F y 
DEN 3 


DEN 3 = iAmI+B") - Q; 


I _ = I- Sin 9 Sin ^ 
yO 0 s s 


I - = I e 

y 1 yO 


'Vl 


q = p, + B‘ + F 
2 y y y 


(28) 

(29) 

(30) 

(31) 

(32) 


P = Jm-' ( p'+2B " ) 
pVb"+p 

L 2 B" 


(33) 

(34) 


M, 




B" 


(35) 


(p,+B"-Q 2 )B + B"F 

DEN 3 


(36) 


i_ y (y=°) = l 2 a 12 + m 2 a 22 + n 2 


(37) 


If the row structure does not exist and plants merge 
into each other y^-*®, A^ 2 =0 an< * A +y (y=y j=®) =0 . However, 

I (y=0) is not necessarily equal to zero. There is no 

transmittance at the end of the row. The reflectance at 
x=0 is given by 



(38) 


I _ x (x=0) = Mj A 21 + N. 


where 


Mj = from axial symmetry 


N 1 " N 2 


A 21 “ P 2 I x0 


T x0 = I 0 Sin 9 s C ° S ♦« 


The input into the sensor depends upon its location rel- 
ative to the position and orientation of the rows. If 
6 are the angular coordinates of the detector then 

' R D ' the input into the detector is given by 

R_ = I , .. S in0_Cos\|f + I, . N S in9 Cost|( 

d -x (x=o) d t d +y(y=y 1 ) d d 


+ I . . Cos 0_ 

- z ( z =0 ) D 


if 2 tt > t D > ~ / 


R = A . .. Sin9 Cosf, + 1 . n .Sin6 Sin if 

D - x ( x = 0 ) D D -y (y =0) D T D 


+ I / Cos 

- z ( z = 0 ) D 


if ~ 7 0; 


R— = A ..Sin9 Sint + I , -.CosQ 
D -y(y=0) D D - z ( z =0 ) D 


if 17 > + I ; 


(39) 

(40) 

(41) 

(42) 


■ (43) 


(44) 

(45) 



+y (y=y 


1 


Sin6 Sin* + A CosO^ . 

D D - z ( z =0) D 


(46) 


if * >TT 
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5 . Shadowing of Soli Between Plants 

In the notation of the last section the plan.t height 

is z. and the solar zenith-angle is 8 . The length of the 
1 s 

shadow is z,tan8 . If it =0 this shadow does not cover the 
1 s s 

open soil between the rows. If iJr^O the length of the 

shadow between the rows is |z 1 tanG S ini|; | . If the total area of 

i s s 

a sensors view is A and — ~ — A, i.e. t 7« of A, is plant cover 

then the solar azimuthal angle^j^O gives rise to an additional 
ground cover, reducing the magnitude of the open soil area tot 


a - a/i t Iz, tan 8 Sin f I -t . 

‘•oil ■ A(1 -Too - -Li — noi LL- > 

100y x 

Therefore, as long as, A >0 the soil reflectance has 

soil 

to be taken into account separately provided also 

| z x tan 0 D Sin * D | < - l)y r 

This later inequality comes from the condition that the 
detector's viewing angles are such that some of the open 
soil is directly 'visible' to the detector. The contri- 
bution of open soil reflectance can be taken into account 


(47) 


by a simple modification of the expression for I 


-z <?z = 0 ) * 


The new expression for I . is denoted by I and is 

- z ( z = 0 ) -z 

given by : 


X -z = (1 


A' (soil ) ) x (R.H.S. of 13) + A ' (soil) x Rg (48) 

A A 


where 

A ' (soil) 


A ( 1 - - )V a n9 s Sint s l 4- \ z 1 tan^S in» p [ ) (49 ) 

100y l/ t 



if.the detector and the sun ace On opposite sides of each other 
i.e. if ^ is in the range of 0 < jji ^ < tt then g is in the 

range g < 2rr . However, if the detector and the sun are such 

that | and ij; g are in the same range of angles 0-tt and rr-2rr 

then I is modified as follows: 


I = (I" ~ )y(R. H .S. Of 13) + Alijoill XRg 


where 


Z 2 


A ' (soil) A ( 1 10Q - 10Q y ) 


(50) 


and is larger of the quantities z ^ tan6 ini|r and 

and z ^ tan 9^ Sin i|r . Implicit in the above expressions 
is the assumption that the soil is a Lambertian reflector. 

For an accurate treatment of this problem one should 
follow Suits and compute contribution from each element of 
the canopy and integrate over the height and width of the 
tree. However, due to lack of time this task will not be 


performed in this report. 
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6 . Critique 

The present model does not contain at least two features that 
should be included to make better predictions from reflec- 
tance. These are: 

(1) Surface reflections from leaves which are smooth 
and may not be Lambertian. 

(2) In the area of a field covered by a plant, there 
are, in general, holes through which sunlight 
falls unattenuated on the ground. At these places 
there is increased ground reflectance and the boun- 
dary condition (10) needs to be modified. 



7 . Experiments that Should Be Performed 

Some of the assumptions made in this theoretical work 
are based on everyday common sense and not on experience 
gained from experiments with plant canopies. In order to 
evaluate the effect of these assumptions or to make as- 
sumptions that are based on true experience we propose 
that experiments should be performed to test the following 
properties and effects in plant canopies: 

(1) Does a plant canopy or a portion of it act as a 
Fresnel reflector? 

(2) Are the scattering and absorption coefficients 

of diffuse light different from those of specular 
light? Do these coefficients depend linearly or 
no n-1 i nearly on the light intensity? 
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